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The asymmetry in the production cross-section a of D mesons, 

_ a(D+)-a(D-) 
P a(D+) + <j(D-y 

is measured in bins of pseudorapidity i] and transverse momentum px within the acceptance 
of the LHCb detector. The result is obtained with a sample of D + — > K®tc + decays 
corresponding to an integrated luminosity of l.Ofb -1 , collected in pp collisions at a 
centre of mass energy of 7 TeV at the Large Hadron Collider. When integrated over the 
kinematic range 2.0 < p^ < 18.0 GeV/c and 2.20 < rj < 4.75, the production asymmetry is 
Ap = (—0.96 ± 0.26 ± 0.18)%. The uncertainties quoted are statistical and systematic, 
respectively. The result assumes that any direct CP violation in the D + — > K^it^ decay is 
negligible. No significant dependence on 77 or p^ is observed. 
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1 Introduction 



The Large Hadron Collider (LHC) offers an excellent opportunity to study heavy flavour 
physics. The rate of production of cc and bb pairs is substantial in the forward region 
close to the beam direction. The associated cross-sections were measured at the LHCb 
experiment in the forward region to be a C c = 1230 ± 190 |J.b and = 74 ± 14 [ib at 
V£=7TeV mi2]. 

Direct production of cc pairs at the LHC occurs almost entirely via QCD and electroweak 
processes that do not discriminate between c and c quarks. However, in hadronization the 
symmetry is broken by the presence of valence quarks, which introduce several processes 
that distinguish between c and c quarks [31 IH E] • For example, a c quark could couple 
to valence quarks to form a charmed baryon, leaving an excess of c quarks. These would 
hadronize to create an excess of D~ mesons over D + mesons. Furthermore, the kinematic 
distributions of charmed hadrons and their antiparticles can differ, introducing production 
asymmetries in local kinematic regions. Analogous production asymmetries in the strange 
sector are well-established at the LHC, and are seen to be large at high rapidity [BJ. 
However, no evidence for a production asymmetry was found in a recent study [7]. 

Searches for CP violation (CPV) in charmed hadron decays can be used to probe 
for evidence of physics beyond the Standard Model [8]. Direct CPV is measured using 
time-integrated observables, and is of particular interest following evidence for CPV in 
two-body D° decays reported by LHCb [9] and subsequently by CDF [ID] . In order to 
understand the origin of this effect, more precise measurements of CP asymmetries in a 
suite of decay modes are required. Production asymmetries have the same experimental 
signature as direct CPV effects and are potentially much larger than the CP asymmetries 
to be determined. This problem can sometimes be avoided by taking the difference in 
asymmetry between two decay modes with a common production asymmetry |9] or by 
studying the difference in kinematic distributions of multi-body decays [TTj. However, 
these methods result in a reduction in statistical power and are not applicable to all final 
states. It is therefore important to measure production asymmetries directly. 

In this Letter, the D ± production asymmetry, defined as 

for cross sections a(D ± ), is determined with a sample of D + — > K®7r + , K® — > 7r + 7r~ 
decays]^] As there are no charged kaons in the final state, the detector biases in this decay 
are simpler to understand than those in other D + decays with higher branching fractions. 
The , a pseudoscalar particle, has a charge-symmetric decay, and the charge asymmetry 
in the pion efficiency at LHCb has been measured previously for the 2011 data sample [7J. 
However, there is the possibility of CPV in the decay. The expected CPV in the D + 
decay, due to the interference of the Cabibbo-favoured and doubly Cabibbo-suppressed 



1 Charge conjugate decays are implied throughout this letter unless stated otherwise. 
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amplitudes, is defined by the charge asymmetry in the partial widths r(D ± ), 

T(D+)-T(D-) 
CP T{D+) + T{D-y 1 J 

Acp is negligible in the Standard Model: a simple consideration of the CKM matrix leads 
to a value of at most 1 x 10~ 4 depending on the strong phase difference between the two 
amplitudes p2]. Since both amplitudes are at tree level, no enhancement of CPV due 
to new physics is expected. The current world-best measurement of Acp, by the Belle 
Collaboration, is consistent with zero: (0.024 ± 0.094 ± 0.067)% [13]. On the other hand, 
CPV in the neutral kaon system induces an asymmetry which must be considered. This 
will be discussed further in Sect. [5] 



2 Detector description 

The LHCb detector [14J is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r] < 5, designed for the study of particles containing b or c quarks. The 
detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector (VELO) surrounding the pp interaction region, a large-area silicon-strip detector 
located upstream of a dipole magnet of reversible polarity with a bending power of 
about 4Tm, and three stations of silicon-strip detectors and straw drift-tubes placed 
downstream. The combined tracking system has a momentum resolution Ap/p that varies 
from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and an impact parameter (IP) resolution 
of 20 |J.m for tracks with high transverse momentum px- Charged hadrons are identified 
using two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and pre-shower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. The 
trigger consists of a hardware stage, based on information from the calorimeter and muon 
systems, an inclusive software stage, which uses the tracking system, and a second software 
stage that exploits the full event information. 



3 Dataset and selection 

The data sample used in this analysis corresponds to l.Ofb -1 of pp collisions taken at a 
centre of mass energy of 7TeV at the Large Hadron Collider in 2011. The polarity of the 
LHCb magnetic field was changed several times during the run, and approximately half 
of the data were taken with each polarity, referred to as 'magnet-up' and 'magnet-down' 
data hereafter. To optimise the event selection and estimate efficiencies, 12.5 million 
pp collision events containing D + — > K®tt + , K® — > tt~tt + decays were simulated with 
Pythia 6.4 [T5] with a specific LHCb configuration [IB]. Decays of hadronic particles are 
described by EvtGen [TTj. The interactions of the generated particles with the detector 
and its response are implemented using the Geant4 toolkit [18] as described in Ref. |19j . 
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Pairs of oppositely charged tracks with a pion mass hypothesis are combined to form 
K® candidates. Only those K® candidates with p^ > 700MeV/c and invariant mass within 
35MeV/c 2 of the nominal value [20j are retained. Surviving candidates are then combined 
with a third charged track, the bachelor pion, to form a D + candidate, with the mass of 
the Kg candidate constrained to its nominal value in a kinematic fit. Each of the three 
pion tracks must be detected in the VELO, so only those K® mesons that decay well 
within the VELO are used. This creates a bias towards short K® decay times. Both the 
Kg and D + candidates are required to have acceptable vertex fit quality. 

Further requirements are applied in order to reduce the background and to align the 
selection of bachelor pions with the dataset used to determine the charge asymmetry in 
the tracking efficiency (see Sect. [6]). The daughters of the K® must have p > 2 GeV/c and 
Pt > 250 MeV/ c. Impact parameter requirements are used to ensure that both the K® 
candidate and its daughter tracks do not originate at any primary vertex (PV) in the event, 
and the K® decay vertex must be at least 10 mm downstream of the PV with which it is 
associated. The bachelor pion must have p > 5 GeV/c and pr > 500 MeV/c, be positively 
identified as a pion rather than as a kaon, electron or muon, and must not come from any 
PV. In addition, fiducial requirements are applied as in Ref. jH] to exclude regions with 
large tracking efficiency asymmetry. All three tracks must have an acceptable track fit 
quality. The D + candidate is required to have pt > 1 GeV/c, to point to a PV (suppressing 
D from B decays), and to have a decay time significantly greater than zero. After these 
criteria are applied, the remaining background is mostly from random combinations of 
tracks. The invariant mass distribution of selected candidates is shown in Fig. [TJ 

In selected events, a trigger decision may be based on part or all of the D + signal 
candidate, on other particles in the event, or both. The second stage of the software 
trigger is required to find a fully reconstructed candidate which meets the criteria to be 
a signal D + — > K®tt + decay. To control potential charge asymmetries introduced by the 
hardware trigger, two possibilities, not mutually exclusive, are allowed. The hardware 
trigger decision must be based on one or both of the K® daughter tracks, or on a particle 
other than the decay products of the D + candidate. In both cases, the inclusive software 
trigger must make a decision based on one of the three tracks that form the D + . For 
the first case, it is explicitly required that the same track activated the hardware trigger, 
and therefore this is independent of the D + charge. The second possibility does not 
depend directly on the D + charge, but an indirect dependence could be introduced if 
the probability for particles produced in association with the signal candidate to activate 
the trigger differs between D + and D~ . This will be discussed further in Sect. [7| After 
applying the selection and trigger requirements, 1,031,068 Kg7i + candidates remain. 

4 Yield determination 

The signal yields are measured in 48 bins of and r\ using binned likelihood fits to 
the distribution of the K®7r + mass m. The bins are shown in Fig. [2j The shapes of the 
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Figure 1: Mass distribution of selected K®7r + candidates. The data are represented by 
symbols with error bars. The dashed curves indicate the signal and the — > K®tt + 
decays, the lower solid line represents the background shape, and the upper solid line 
shows the sum of all fit components. 

Df g -. — > K®ti + mass peaks are described by 'CruijfF functions [2"T] . 



with the measured masses defined by the free parameter /z, the widths by <jl and or, and 
the tails by and a^. The parameters and are used for m < fi and ccr and a R 
for m > fi. The background is fitted with a straight line plus an additional Gaussian 
component to account for background from — > JC?7r + 7r° decays. The yield of the latter 
is consistent with zero in most px, V bins. The fit is performed simultaneously over four 
subsamples (D + magnet-up, D + magnet-down, D~ magnet-up, and D~ magnet-down 
data) with the masses and yields of the D^, and the yield of background, allowed to 
vary independently in the four subsamples. All other parameters are shared. The charge 
asymmetries are then determined from the yields. The results are cross-checked with a 
sideband subtraction procedure under the assumption of a linear background. 

5 Effect of CP violation in the neutral kaon system 

CP violation in the neutral kaon system can affect the observed asymmetry in the 
D + — > Kgii + decay [22]. The bias on Ap due to the CPV depends on the decay time 
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Figure 2: Background-subtracted distribution of transverse momenta px versus pseudora- 
pidity 7] for selected D + — > K^tt + candidates in a signal region of 1845 < m < 1890 MeV/c 2 . 
The bin marked with an asterisk is excluded from the weighted average over the production 
asymmetries in the bins used to obtain the final result. 



acceptance F(t) of the K® meson, according to 

J °° F(t)e-^ r s+^ (cos Amt - |g sin Ami) dt 



A e ~ 23?(e) 



- r stdt 



(4) 



where e parameterises the indirect CPV in neutral kaon mixing, Tg and IY are the decay 
widths of the K® and K® respectively, and Am is their mass difference [23l l2l] . Direct 
CPV and terms of order e 2 are neglected. To determine the decay time acceptance, 
the decay time is fitted with an empirical function shown in Fig. |3j All of the K® 
candidates used in this analysis decay inside the VELO with an average measured lifetime 



of 6.97 ± 0.02 ps, which is much shorter than the nominal K~ lifetime of 89.5 ps. Using 
E(e) = 1.65 x 10" 3 [20J in Eq. g we obtain A e = (2.831^;°°^ 



1.004^ 



x 10^ 4 for the CPV in 

the neutral kaon system, where the uncertainty quoted is statistical only. This value is 
subtracted from the measured production asymmetry and a systematic uncertainty equal 
to its central value is assigned. 



6 Results 

In order to convert the measured charge asymmetries in the 48 bins of pr and r\ into 
production asymmetries, a correction for the asymmetry in the pion reconstruction efficiency 
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Figure 3: Observed K® decay time distribution within the LHCb acceptance. The data 
points are fitted with an empirical function (solid curve). This contains a component for 
the upper decay time acceptance, due mainly to the requirement that the K® decays inside 
the VELO (dashed curve) and a component for the lower decay time acceptance, due to 
the selection cuts (dotted curve). These are shown scaled by arbitrary factors. The CPV 
is not sensitive to the fine details of the distribution, so the fit quality is not important. 



is made. This asymmetry was evaluated previously in eight bins of pion azimuthal angle 
and two bins of pion momentum with a control sample of D* + — > D tt + , D° — > K~7t + 7t~tt + 
decays in the same dataset [7j. The average efficiency asymmetry ratios e^+fe^- in that 
sample were found to be 0.9914 ± 0.0040 for magnet-up data and 1.0045 ± 0.0034 for 
magnet-down data. 

After the correction is applied, the resulting asymmetries for magnet-up and magnet- 
down data in each D + px and rj bin are averaged with equal weights to obtain the 
production asymmetries in two-dimensional bins of px and r/, given in Table[TJ Any left-right 
asymmetries that differ between the signal D + — > K®tt + decay and the D° — > K~7t + 7t~tt + 
control channel will cancel in this average. 

Reconstruction and selection efficiencies from the simulation are used to calculate 
binned efficiency-corrected yields. These are used to weight the production asymmetries 
in the average over the px and i] bins. The result is an asymmetry for D + produced in 
the LHCb acceptance. The same weighting technique is applied to obtain production 
asymmetries as one-dimensional functions of px and 77. The bin marked with an asterisk 
in Fig. [2] has a high cross section but is mostly outside the acceptance and so it is excluded 
from the average. After subtracting the contribution from CPV in the kaon system, the 
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Figure 4: Production asymmetry as a function of (a) transverse momentum p^ and (b) 
pseudorapidity rj. The straight line fits have slopes of (0.09 ± 0.07) x 10~ 2 ( GeV/c) -1 and 
(—0.36 ± 0.28)%, and values of \ 2 P er degree of freedom of 5.5/6 and 2.2/4, respectively. 
The error bars include only the statistical uncertainty on the D + signal sample and are 
uncorrelated within a given plot. 



Table 1: Production asymmetry for D + mesons, in percent, in (pt,t)) bins, for 2.0 < 
Pt < 18.0GeV/c and 2.20 < r\ < 4.75. The uncertainties shown are statistical only; the 
systematic uncertainty is 0.17% (see Table [2]). 



p T ( GeV/c) 


(2.20,2.80) 


(2.80,3.00) 


rj 

(3.00,3.25) 


(3.25,3.50) 


(3.50,3.80) 


(3.80,4.75) 


(2.00,3.20) 


-0.0 ±2.5 


-2.2 ±1 


2 


-0.4 ±0.8 


-0.4 ±0.7 


-1.2 ±0.6 


-1.2 ±0.5 


(3.20,4.00) 


-0.4 ±0.9 


-0.4 ±0 


7 


-0.4 ±0.5 


-1.1 ±0.5 


±0.1 ±0.5 


-1.2 ±0.5 


(4.00,4.55) 


±0.1 ±0.8 


-1.0 ±0 


8 


-1.3 ±0.6 


-2.0 ±0.6 


-0.1 ±0.6 


-2.1 ±0.7 


(4.55,5.20) 


-1.6 ±0.7 


-0.6 ±0 


8 


-0.5 ±0.6 


-0.7 ±0.6 


-1.6 ±0.6 


-2.0 ±0.8 


(5.20,6.00) 


-0.5 ±0.7 


-0.8 ±0 


8 


±0.2 ±0.7 


-0.3 ±0.7 


-0.6 ±0.7 


-1.2 ±0.9 


(6.00,7.00) 


-1.4 ±0.8 


±0.5 ±1 





-0.9 ±0.9 


-0.6 ±0.9 


-0.7 ±0.9 


-1.6 ±1.2 


(7.00,9.50) 


-0.4 ±0.8 


-0.4 ± 1 


1 


-0.2 ±1.1 


±1.7 ± 1.1 


-1.4 ± 1.1 


±1.2 ±1.4 


(9.50, 18.00) 


-0.6 ±1.3 


±1.8 ±2 


3 


-2.5 ±2.2 


±1.8 ±2.4 


±1.1 ±2.5 


-7± 11 



production asymmetry is (—0.96 ± 0.19 ± 0.18)%. The uncertainties are the statistical 
errors on the D + — > K®7r + yields and that due to the tagged D° — > K~7t + 7t~tt + sample 
used to calculate the pion efficiencies. Summing these in quadrature, we obtain 

A P = (-0.96 ± 0.26 (stat.)) %. 

The production asymmetry as a function of pr and rj is given in Fig. |4j No significant 
dependence of the asymmetry on these variables is observed. As a cross-check, the average 
production asymmetry is calculated for magnet-up and magnet-down data separately, and 
found to be fully consistent: (-1.07 ± 0.41)% and (-0.85 ± 0.34)%, respectively. 
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Table 2: Summary of absolute values of systematic uncertainties on Ap. For the binned 
production asymmetries given in Table [TJ all uncertainties except that on the reconstruction 
efficiency apply, giving a combined systematic uncertainty of 0.17%. 



Systematic effect Uncertainty (%) 

Trigger asymmetries 0.15 

D from B 0.04 

Selection criteria 0.05 

Running conditions 0.04 

Pion efficiency 0.02 

Fitting 0.04 

Kaon CP violation 0.03 

Weights (reconstruction efficiency) 0.05 

Total including uncertainty on weights 0.18 



7 Systematic uncertainties 

The sources of systematic uncertainty are summarised in Table [2j The dominant uncertainty 
of 1.5 x 10~ 3 is due to asymmetries introduced by the trigger. Events which are triggered 
independently of the signal decay, i.e. by a track that does not form part of the signal 
candidate, could be triggered by particles produced in association with the D + meson. 
If this occurs, the asymmetry in this sample would be correlated with the production 
asymmetry, and would bias the measurement of it. This was studied with a control 
sample of the abundant D + — > K~tt + tt + decay. To mimic the charge-unbiased sample of 
D + —7- Kg7r + decays which are triggered by a K® daughter, we choose the kaon and one 
pion at random and require that the trigger decision is based on one of these tracks. This 
is close to being charge-symmetric between D + and D~ candidates, with some residual 
effects due to differences in material interaction between K + and K~ mesons. The raw 
asymmetry in this subsample of D + — > K~7t + tt + decays is then compared to that in the 
much larger sample of candidates that are triggered independently of the signal decay. 
The difference in raw charge asymmetry between these two samples, (1.5 ± 0.4) x 10 -3 , is 
a measure of the scale of the bias. Unlike the signal, the K~tt + tt + decay also includes a 
component due to the K + /K~ asymmetry, and therefore this is treated as a systematic 
uncertainty rather than a correction. This is cross checked with other control samples 
such as -D+ — > 07r + and the uncertainty is found to be conservative. 

Further systematic uncertainties arise from the contamination of the prompt sample 
by D candidates that originate from B decays. The yield of these is calculated using 
the measured cross-sections [TJ [2] , branching ratios, and efficiencies determined from the 
simulation. The fraction of D candidates from B decays is found to be (1.2 ± 0.3)%. 
This quantity is combined with the B° production asymmetry, which is estimated to be 
(—1.0 ± 1.3)% [25J, to determine the systematic uncertainty. 



S 



Certain selection criteria differ between the D + — > K®tt + signal sample and the 
D° K~7t + 7t~7t + decays used to determine the asymmetry in the pion efficiencies. The 
charge asymmetry is found to depend weakly on the value of the requirement on the pion 
Pt- Pions in the signal sample must have pt > 500 MeV/c while those in the control sample 
must have pt > 300MeV/c. A systematic uncertainty is calculated by estimating the 
proportion of signal candidates with 300 < < 500 MeV/ c and multiplying this fraction 
by the difference between the charge asymmetries in the low px region and the average. 

The difference in signal yields per pb -1 of integrated luminosity between magnet-up 
and magnet-down data is used to determine a systematic uncertainty for changes in 
running conditions that could impair the cancellation of detector asymmetries achieved by 
averaging over the magnet polarities. There is also a systematic uncertainty on the pion 
efficiency asymmetry associated with the determination of the yields of D° — > K~7t + 7t~tt + 
decays. The error associated with the mass fit is determined by comparing fitted and 
sideband-subtracted results. The CPV in the neutral kaon decay, discussed in Sect. |5j is 
also included as a systematic uncertainty. 

Other systematic effects such as regeneration in the neutral kaon system [26] . second 
order effects due to the kinematic binning of the D + — > K®7t + sample, and asymmetric 
backgrounds such as that from — > K®K + with the kaon misidentified as a pion, were 
considered but found to be negligible. When taking the average asymmetry weighted by 
the efficiency-corrected yield in each bin, the limited number of simulated events leads to 
an uncertainty on the reconstruction efficiency and hence on the per-bin weights. This 
does not contribute to the uncertainty on the individual asymmetries given in Table [TJ 
which are calculated without using the simulation. A quadratic sum yields an overall 
systematic uncertainty of 1.8 x 10~ 3 . 

In principle, CPV in the charm decay could occur via the interference of Cabibbo- 
favoured and doubly Cabibbo-suppressed amplitudes, but this is strongly suppressed by 
the CKM matrix and no evidence for it has been observed at the B- factories [27] [13] . If 
we allowed for the possibility of new physics or large unexpected enhancements of the 
Standard Model CPV in these tree-level D + decays, the uncertainty on the null result 
found at Belle [13] would increase the total systematic uncertainty to 2.1 x 10 -3 . 

8 Conclusions 

Evidence for a charge asymmetry in the production of D + decays is observed at LHCb. 
In the kinematic range 2.0 < px < 18.0GeV/c and 2.20 < rj < 4.75, excluding the region 
with 2.0 < pt < 3.2 GeV/c, 2.20 < r] < 2.80, the average asymmetry is 

A P = (-0.96 ± 0.26 ± 0.18)%, 

where the first uncertainty is statistical and the second is systematic. The result is 
inconsistent with zero at approximately three standard deviations. There is no evidence 
for a significant dependence on p^< or pseudorapidity at the present level of precision. 
The bias on the measured asymmetry due to CP violation in kaon decays has been 
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calculated and found to be almost negligible for this dataset. These results are consistent 
with expectations [5] and lay the foundations for searches for CP violation in Cabibbo 
suppressed D + decays. 
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